Deep population structure linked to host vernalization requirement in the barley net blotch fungal pathogen.

Invasive fungal pathogens pose a substantial threat to widely cultivated crop species, owing to their capacity to adapt to new hosts and new environmental conditions. Gaining insights into the demographic history of these pathogens and unravelling the mechanisms driving coevolutionary processes are crucial for developing durably effective disease management programmes. Pyrenophora teres is a significant fungal pathogen of barley, consisting of two lineages, Ptt and Ptm, with global distributions and demographic histories reflecting barley domestication and spread. However, the factors influencing the population structure of P. teres remain poorly understood, despite the varietal and environmental heterogeneity of barley agrosystems. Here, we report on the population genomic structure of P. teres in France and globally. We used genotyping-by-sequencing to show that Ptt and Ptm can coexist in the same area in France, with Ptt predominating. Furthermore, we showed that differences in the vernalization requirement of barley varieties were associated with population differentiation within Ptt in France and at a global scale, with one population cluster found on spring barley and another population cluster found on winter barley. Our results demonstrate how cultivation conditions, possibly associated with genetic differences between host populations, can be associated with the maintenance of divergent invasive pathogen populations coexisting over large geographic areas. This study not only advances our understanding of the coevolutionary dynamics of the Pt-barley pathosystem but also prompts further research on the relative contributions of adaptation to the host versus adaptation to abiotic conditions in shaping Ptt populations.

Quantitative and qualitative plant-pathogen interactions call upon similar pathogenicity genes with a spectrum of effects.

Septoria leaf blotch is a foliar wheat disease controlled by a combination of plant genetic resistances and fungicides use. R-gene-based qualitative resistance durability is limited due to gene-for-gene interactions with fungal avirulence (Avr) genes. Quantitative resistance is considered more durable but the mechanisms involved are not well documented. We hypothesize that genes involved in quantitative and qualitative plant-pathogen interactions are similar. A bi-parental population of Zymoseptoria tritici was inoculated on wheat cultivar 'Renan' and a linkage analysis performed to map QTL. Three pathogenicity QTL, Qzt-I05-1, Qzt-I05-6 and Qzt-I07-13, were mapped on chromosomes 1, 6 and 13 in Z. tritici, and a candidate pathogenicity gene on chromosome 6 was selected based on its effector-like characteristics. The candidate gene was cloned by Agrobacterium tumefaciens-mediated transformation, and a pathology test assessed the effect of the mutant strains on 'Renan'. This gene was demonstrated to be involved in quantitative pathogenicity. By cloning a newly annotated quantitative-effect gene in Z. tritici that is effector-like, we demonstrated that genes underlying pathogenicity QTL can be similar to Avr genes. This opens up the previously probed possibility that 'gene-for-gene' underlies not only qualitative but also quantitative plant-pathogen interactions in this pathosystem.

Complete Genome Sequences of Septoria linicola: A Resource for Studying a Damaging Flax Pathogen.

Fungal genus Septoria causes diseases in a wide range of plants. Here, we report the first genome sequences of two strains of Septoria linicola, the causal agent of the pasmo disease of flax (Linum usitatissimum). The genome of the first strain, SE15195, was fully assembled in 16 chromosomes, while 35 unitigs were obtained for a second strain, SE14017. Structural annotations predicted 13,096 and 13,085 protein-encoding genes and transposable elements content of 19.0 and 18.1% of the genome for SE15195 and SE14017, respectively. The four smaller chromosomes 13 to 16 show genomics features of potential accessory chromosomes. The assembly of these two genomes is a new resource for studying S. linicola and improving management of pasmo. [Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

First report of barley virus G infecting winter barley (Hordeum vulgare L.) in France.

As part of a cereals virome project high throughput sequencing (HTS)-based viral indexing was performed on plants with symptoms of barley yellow dwarf disease collected in June (2017-2020) in the main French cereals production areas. Total RNAs from 32 individual plants were purified (RNeasy Plant Mini Kit, Qiagen, Courtaboeuf, France) and Illumina sequenced (2x150 nt) following ribodepletion (Genewiz-Azenta, Leipzig, Germany). Following quality trimming, reads for each sample were de novo assembled (CLC Genomics Workbench 21, Qiagen) [1] and contigs annotated by BlastX analysis. In four winter barley samples collected in 2018 (18-58, 18-325 and 18-326) and 2019 (19-30A), besides contigs representing diverse viruses such as barley yellow dwarf viruses-PAV and PAS, Hordeum vulgare endornavirus, cereal yellow dwarf virus-RPV (18-326), wheat dwarf virus (18-325 and 18-326) and a novel Polerovirus (18-58 and 18-326), large contigs with high identity to barley virus G (BVG) were identified. BVG, a tentative Polerovirus, was initially reported in barley in South Korea in 2016 [2] and has so far been identified in a few other hosts including wheat, oat, maize, proso and foxtail millets as well as switchgrass. It has been reported from the USA and Australia [3] and, in Europe, from the Netherlands, Germany, Hungary and Greece [4]. Large BVG scaffolds representing near complete genomes could be reconstructed for each sample, integrating a total of 128.339, 7.188, 8.078 and 20.073 reads, for samples 19-30A, 18-58, 18325 and 18-326 respectively. Given that between 17.2 and 20.5 million reads had been obtained per sample, these values translate into between 0.04% (18-58 and 18-325) and 0.6% (19-30A) of total reads, and to average coverages of between 158x (18-58) and 2866x (19-30A) for the genomic scaffolds. The four assembled sequences (5584-5610 nt) have been deposited in GenBank (ON419453-ON419456). They are nearly identical (98.4 to 99.5% nt identity) and share between 97.7% and 98.5% nt identity with a barley reference isolate from the South Korea (NC_029906). To confirm the presence of BVG, a primer pair was designed based on available BVG sequences. Primers BVG-F(5'-CTAGCCCAACGAGTTGCGGG-3') and BVG-R(5'-GGTACAGAAGCTCTACGGTTC-3') amplifying a 394 nt were used in a two-step RT-PCR on new RNA extracts obtained from the 18-325 and 18-326 infected plants. The amplicons were directly sequenced and showed respectively 99.2% (ON419457, 18-325) and 100% (18-326) nt identity with the corresponding de novo scaffolds. The four analyzed samples have been collected respectively in 2018 (18-58, 18-325, 18-326) and 2019 (19-30A) in three different regions of France (Auvergne-Rhône-Alpes, Occitanie and Centre-Val de Loire), indicating a wide distribution and a persistence over time of BVG in France. To our knowledge, this represents the first report of a natural infection of BVG in cultivated winter barley in France. Presence of BVG may have been overlooked in a range of situation, as indicated by its retrospective discovery in a 34 years old Australian sample [3], possibly explaining its broad distribution in France. While the mixed infection status of the analyzed plants precludes any conclusion on its pathogenicity in French cereals, BVG has been reported to be associated with a range of symptoms in various hosts so that further studies to evaluate its prevalence and impact in France and to begin to understand its epidemiology are clearly warranted by the present results.

Wheat Stem Rust Back in Europe: Diversity, Prevalence and Impact on Host Resistance.

The objective of this study was to investigate the re-emergence of a previously important crop pathogen in Europe, Puccinia graminis f.sp. tritici, causing wheat stem rust. The pathogen has been insignificant in Europe for more than 60 years, but since 2016 it has caused epidemics on both durum wheat and bread wheat in local areas in southern Europe, and additional outbreaks in Central- and West Europe. The prevalence of three distinct genotypes/races in many areas, Clade III-B (TTRTF), Clade IV-B (TKTTF) and Clade IV-F (TKKTF), suggested clonal reproduction and evolution by mutation within these. None of these genetic groups and races, which likely originated from exotic incursions, were detected in Europe prior to 2016. A fourth genetic group, Clade VIII, detected in Germany (2013), was observed in several years in Central- and East Europe. Tests of representative European wheat varieties with prevalent races revealed high level of susceptibility. In contrast, high diversity with respect to virulence and Simple Sequence Repeat (SSR) markers were detected in local populations on cereals and grasses in proximity to Berberis species in Spain and Sweden, indicating that the alternate host may return as functional component of the epidemiology of wheat stem rust in Europe. A geographically distant population from Omsk and Novosibirsk in western Siberia (Russia) also revealed high genetic diversity, but clearly different from current European populations. The presence of Sr31-virulence in multiple and highly diverse races in local populations in Spain and Siberia stress that virulence may emerge independently when large geographical areas and time spans are considered and that Sr31-virulence is not unique to Ug99. All isolates of the Spanish populations, collected from wheat, rye and grass species, were succesfully recovered on wheat, which underline the plasticity of host barriers within P. graminis. The study demonstrated successful alignment of two genotyping approaches and race phenotyping methodologies employed by different laboratories, which also allowed us to line up with previous European and international studies of wheat stem rust. Our results suggest new initiatives within disease surveillance, epidemiological research and resistance breeding to meet current and future challenges by wheat stem rust in Europe and beyond.

Resistance of the Wheat Cultivar 'Renan' to Septoria Leaf Blotch Explained by a Combination of Strain Specific and Strain Non-Specific QTL Mapped on an Ultra-Dense Genetic Map.

Quantitative resistance is considered more durable than qualitative resistance as it does not involve major resistance genes that can be easily overcome by pathogen populations, but rather a combination of genes with a lower individual effect. This durability means that quantitative resistance could be an interesting tool for breeding crops that would not systematically require phytosanitary products. Quantitative resistance has yet to reveal all of its intricacies. Here, we delve into the case of the wheat/Septoria tritici blotch (STB) pathosystem. Using a population resulting from a cross between French cultivar Renan, generally resistant to STB, and Chinese Spring, a cultivar susceptible to the disease, we built an ultra-dense genetic map that carries 148,820 single nucleotide polymorphism (SNP) markers. Phenotyping the interaction was done with two different Zymoseptoria tritici strains with contrasted pathogenicities on Renan. A linkage analysis led to the detection of three quantitative trait loci (QTL) related to resistance in Renan. These QTL, on chromosomes 7B, 1D, and 5D, present with an interesting diversity as that on 7B was detected with both fungal strains, while those on 1D and 5D were strain-specific. The resistance on 7B was located in the region of Stb8 and the resistance on 1D colocalized with Stb19. However, the resistance on 5D was new, so further designated Stb20q. Several wall-associated kinases (WAK), nucleotide-binding and leucine-rich repeats (NB-LRR) type, and kinase domain carrying genes were present in the QTL regions, and some of them were expressed during the infection. These results advocate for a role of Stb genes in quantitative resistance and for resistance in the wheat/STB pathosystem being as a whole quantitative and polygenic.

Fusarium species and enniatin mycotoxins in wheat, durum wheat, triticale and barley harvested in France.

Contamination with enniatins A, A1, B and B1 was investigated in 1240 samples of small grain cereals (470 wheat, 260 durum wheat, 282 spring barley, 172 triticale and 56 winter barley) from the French harvests of 2012 to 2014. Associations with Fusarium avenaceum, F. tricinctum and F. poae were assessed, with the identification of Fusarium species by real-time PCR and mycotoxin quantification by LC-MS/MS. Enniatins were common in the fields sampled. Enniatin concentrations varied between years but were consistently highest on spring barley (mean values of 199 to 1316 µg/kg) and triticale (mean values from 131 to 1218 µg/kg), and lower on wheat (mean values from 47 to 142 µg/kg) and durum wheat (mean values from 55 to 596 µg/kg). The concentrations of the various enniatins were strongly correlated with each other (Pearson's correlation coefficient of 0.61 to 0.98). Enniatin B was the most frequent (68% of the total enniatin content), followed by enniatin B1 (22%), enniatin A1 (7%) and enniatin A (3%). Fusarium species were quantified by calculating arithmetic mean total DNA levels. F. tricinctum was the most abundant (0.177 pg/ng total DNA), followed by F. avenaceum (0.141 pg/ng total DNA) and F. poae (0.091 pg/ng total DNA). Total DNA levels for each species, and the predominant species varied between years and crops. Small grain cereal species (p value < 0.001), harvest year (p value = < 0.001) and the presence of F. avenaceum (p value < 0.001), F. tricinctum (p value < 0.001) or F. poae (p value = 0.017) affected enniatin content. F. tricinctum was the leading enniatin producer on durum wheat (29% to 45%) and spring barley (23 to 37%). F. avenaceum produced large amounts of enniatins on durum wheat (13% to 17%) and wheat (1% to 18%) and was the leading producer on triticale (30% to 55%). F. poae made a minor contribution on all crops, never accounting for more than 2% of total enniatin content. Enniatins are, thus, highly prevalent in French small grain cereals and are mostly produced by F. avenaceum and F. tricinctum.

Metabarcoding targeting the EF1 alpha region to assess Fusarium diversity on cereals.

Fusarium head blight (FHB) is a major cereal disease caused by a complex of Fusarium species. These species vary in importance depending on climatic conditions, agronomic factors or host genotype. In addition, Fusarium species can release toxic secondary metabolites. These mycotoxins constitute a significant food safety concern as they have health implications in both humans and animals. The Fusarium species involved in FHB differ in their pathogenicity, ability to produce mycotoxins, and fungicide sensitivity. Accurate and exhaustive identification of Fusarium species in planta is therefore of great importance. In this study, using a new set of primers targeting the EF1α gene, the diversity of Fusarium species on cereals was evaluated using Illumina high-throughput sequencing. The PCR amplification parameters and bioinformatic pipeline were optimized with mock and artificially infected grain communities and further tested on 65 field samples. Fusarium species were retrieved from mock communities and good reproducibility between different runs or PCR cycle numbers was be observed. The method enabled the detection of as few as one single Fusarium-infected grain in 10,000. Up to 17 different Fusarium species were detected in field samples of barley, durum and soft wheat harvested in France. This new set of primers enables the assessment of Fusarium diversity by high-throughput sequencing on cereal samples. It provides a more exhaustive picture of the Fusarium community than the currently used techniques based on isolation or species-specific PCR detection. This new experimental approach may be used to show changes in the composition of the Fusarium complex or to detect the emergence of new Fusarium species as far as the EF1α sequence of these species show a sufficient amount of polymorphism in the portion of sequence analyzed. Information on the distribution and prevalence of the different Fusarium species in a given geographical area, and in response to various environmental factors, is of great interest for managing the disease and predicting mycotoxin contamination risks.

Low Amplitude Boom-and-Bust Cycles Define the Septoria Nodorum Blotch Interaction.

INTRODUCTION: Septoria nodorum blotch (SNB) is a complex fungal disease of wheat caused by the Dothideomycete fungal pathogen Parastagonospora nodorum. The fungus infects through the use of necrotrophic effectors (NEs) that cause necrosis on hosts carrying matching dominant susceptibility genes. The Western Australia (WA) wheatbelt is a SNB "hot spot" and experiences significant under favorable conditions. Consequently, SNB has been a major target for breeders in WA for many years. MATERIALS AND METHODS: In this study, we assembled a panel of 155 WA P. nodorum isolates collected over a 44-year period and compared them to 23 isolates from France and the USA using 28 SSR loci. RESULTS: The WA P. nodorum population was clustered into five groups with contrasting properties. 80% of the studied isolates were assigned to two core groups found throughout the collection location and time. The other three non-core groups that encompassed transient and emergent populations were found in restricted locations and time. Changes in group genotypes occurred during periods that coincided with the mass adoption of a single or a small group of widely planted wheat cultivars. When introduced, these cultivars had high scores for SNB resistance. However, the field resistance of these new cultivars often declined over subsequent seasons prompting their replacement with new, more resistant varieties. Pathogenicity assays showed that newly emerged isolates non-core are more pathogenic than old isolates. It is likely that the non-core groups were repeatedly selected for increased virulence on the contemporary popular cultivars. DISCUSSION: The low level of genetic diversity within the non-core groups, difference in virulence, low abundance, and restriction to limited locations suggest that these populations more vulnerable to a population crash when the cultivar was replaced by one that was genetically different and more resistant. We characterize the observed pattern as a low-amplitude boom-and-bust cycle in contrast with the classical high amplitude boom-and-bust cycles seen for biotrophic pathogens where the contrast between resistance and susceptibility is typically much greater. Implications of the results are discussed relating to breeding strategies for more sustainable SNB resistance and more generally for pathogens with NEs.

Classical and next generation sequencing approaches unravel Bymovirus diversity in barley crops in France.

Despite the generalized use of cultivars carrying the rym4 resistance gene, the impact of viral mosaic diseases on winter barleys increased in recent years in France. This change could reflect i) an increased prevalence of the rym4 resistance-breaking pathotype of Barley yellow mosaic virus Y (BaYMV-2), ii) the emergence of rym4 resistance-breaking pathotypes of Barley mild mosaic virus (BaMMV) or iii) the emergence of other viruses. A study was undertaken to determine the distribution and diversity of viruses causing yellow mosaic disease. A collection of 241 symptomatic leaf samples from susceptible, rym4 and rym5 varieties was gathered from 117 sites. The viruses present in all samples were identified by specific RT-PCR assays and, for selected samples, by double-stranded RNA next-generation sequencing (NGS). The results show that BaYMV-2 is responsible for the symptoms observed in varieties carrying the resistance gene rym4. In susceptible varieties, both BaYMV-1 and BaYMV-2 were detected, together with BaMMV. Phylogenetic analyses indicate that the rym4 resistance-breaking ability independently evolved in multiple genetic backgrounds. Parallel analyses revealed a similar scenario of multiple independent emergence events in BaMMV for rym5 resistance-breaking, likely involving multiple amino acid positions in the viral-linked genome protein. NGS analyses and classical techniques provided highly convergent results, highlighting and validating the power of NGS approaches for diagnostics and viral population characterization.